Trapping ultracold dysprosium: a highly magnetic gas for dipolar physics 
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Ultracold dysprosium gases, with a magnetic moment ten times that of alkali atoms and equal only 
to terbium as the most magnetic atom, are expected to exhibit a multitude of fascinating coUisional 
dynamics and quantum dipolar phases, including quantum liquid crystal physics. We report the 
first laser cooling and trapping of half a billion Dy atoms using a repumper-free magneto-optical 
trap (MOT) and continuously loaded magnetic confinement, and we characterize the trap recycling 
dynamics for bosonic and fermionic isotopes. The first inelastic collision measurements in the few 
partial wave, 100 mK, regime are made in a system possessing a submerged open electronic 

f-shell. In addition, we observe unusual stripes of intra-MOT < 10 /xK sub-Doppler cooled atoms. 

PACS numbers: 37.10.De, 37.10.Gh, 37.10.Vz, 71.10.Ay 



Ultracold gases of extraordinarily magnetic atoms, 
such as dysprosium, offer opportunities to explore 
strongly correlated matter in the presence of the long- 
range, anisotropic dipole-dipole interaction (DDI). Such 
interactions in the presence (or absence) of polarizing 
fields can compete with short-range interactions to in- 
duce phases beyond those described by the nearest neigh- 
bor Hubbard model [T]. Specifically, quantum liquid 
crystal (QLC) physics (see Ref. [2j and citations within) 
describes strongly correlated systems in which a Fermi 
surface can spontaneously distort (nematics) or cleave 
into stripes (smectics) 3 . While material complexity 
can inhibit full exploration of QLC phases in condensed 
matter, QLC phases may be more extensively character- 
ized in tunable ultracold gases. In contrast to ultracold 
ground state polar molecules ultracold Dy offers the 
ability to explore the spontaneously broken symmetries 
inherent in QLCs since the DDI is realized without a 
polarizing field. An exciting prospect lies in observing 
spontaneous magnetization in dipolar systems, e.g., the 
existence of a quantum ferro-nematic phase in ultracold 
fermionic Dy gases not subjected to a polarizing field [5] . 

We report the first magneto-optical trap (MOT) and 
ultracold coUisional rates of this highly complex atom. 
The stable atoms possessing the largest magnetic mo- 
ments are the neighboring lanthanide rare-earths, Tb and 
Dy (both 10 Bohr magnetons {hb) to within 0.6% [6l[7]). 
Prior to the present work, the coldest Dy temperatures 
were achieved via buffer gas and adiabatic cooling to 
50 mK with final densities of < 10^ cm~^ [5], and Dy 
beams have been transversely pushed and unidirection- 
ally cooled via photon scattering [3 [TU] . 

Recent experiments using degenerate ^^Cr, a bosonic 
S-state atom with 6 of magnetic moment, have be- 
gun to explore quantum ferrofluids [T^]. With suitable 
scattering lengths, the larger magnetic moment — and 9x 
larger DDI-mass ratio — of Dy should allow experimental 
access beyond the superfluid and Mott insulator regions 
of the extended Bose-Hubbard phase diagram to the den- 
sity wave and supersolid regimes [T3]. Co-trapping iso- 
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FIG. 1. (color online), (a) Dysprosium energy level struc- 
ture [H] [TT]. The MOT and Zeeman slower employ the 
strongest laser cooling transition between the even parity 
(red) ground state and the odd (black) excited state at 421 
nm. Dy has five high-abundance isotopes; three bosons 
(164 I62j3y^ leopy ^^^^-^ ^ ^ femiious ("^^Dy 

and ^^^Dy with / — 5/2). (b, c) Fermion hyperfine structure 
shown along with the cooling transition and energy splittings 
(not to scale). Poor optical pumping reduces the '^^^Dy MOT 
population due to a large 14r frequency difference between 
the f =21/2^F'=23/2 and F=19/2-»F'=21/2 transitions. 

topes of Dy will allow exploration of dipolar Bose-Fermi 
mixtures of near equal mass. Studies of degenerate spinor 
gases |14j with large spins, and simulations of dense nu- 
clear matter |15j are further exciting avenues of research. 

In addition, ultracold samples of Dy will aid precision 
measurements of parity nonconservation and variation of 
fundamental constants [S], single-ion implantation |16j . 
and quantum information science [17'. The low- lying 
telecommunications band (1322 nm) and InAs quantum 
dot (QD) amenable (1001 nm) transitions (Fig.[T]) will en- 
able hybrid quantum circuits of atom-photonic or atom- 
QD systems. Novel coUisional physics [18 and complex 
molecular association phenomena are expected in clouds 
of these L 7^ atoms possessing submerged open f-shells 
inside closed 5s and 6s electron shells. 
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FIG. 2. (color online), (a) Dy MOT recycling and contin- 
uously loaded magnetic trap (MT) schematic, (b) Popula- 
tion ratio iVivroT of recaptured MOT to steady state MOT. 
Black line is fit using Eqs. [l] to TVmot with Zeeman slower 
and atomic beam off, and tm = 1 s delay between steady 
state MOT and recapture. Inset (c) plots the peak recap- 
tured MOT population versus hold time in MT. Lines are 
fits to t < f s data with exponential time constants fixed at 
-^siow ~ l-"^ ^ and -Rgiow ~ 2.3 s 

Dy structure is similar to the recently magneto- 
optically trapped Er (7 ^s) [IS]- Despite the multitude 
of excited state population loss channels, the repumper- 
less Dy MOT forms in a similar fashion to Er's: The 
large magnetic moment allows Dy to remain magnetically 
trapped in the MOT region while excited state popula- 
tion decays through the metastable states. No repump- 
ing lasers are necessary since a sufficient fraction of the 
atoms recycle to the MOT to overcome loss. 

The Dy cooling and trapping apparatus consists of a 
high temperature oven, transverse cooling stage, Zeeman 
slower, and MOT trapping region. We elaborate on the 
design elsewhere [12] with brief details here. Nuggets 
of Dy are heated to 1250 °C in a Ta crucible aper- 
tured to provide an atomic beam. Three ion pumps — 
located at the oven, transverse pumping stage, and MOT 
chamber — along with titanium sublimation, achieve a 
vacuum of 1.8 X 10~" Torr during MOT operation. 

The atomic beam passes through a differential pump- 
ing tube before intersecting four 100 mW, 2 cm-long 
transverse cooling beams detuned -0.2r from the broad 
421-nm transition [20 . Transverse cooling enhances the 
optimized MOT population by a factor of 3-4. The 
atomic beam then passes through a spin-flip Zeeman 
slower [10' operating at -24r detuning from the 421- 
nm line. The power of the Zeeman slowing laser is >1 
W at the slower entrance; we find a linear increase in 
MOT population until ~1 W. Atoms in the slowed beam 
are captured by a 6-beam MOT with detuning S = — F 
from the 421-nm line and 12 mW in each 1.1-cm e^^- 
waist beam (total intensity 7 = 36 mW/cm^). We ob- 
serve maximum MOT population for a 35 G/cm mag- 
netic quadrupole gradient VzB « 2VxB. The 421-nm 
lasers are derived from two doubled Ti:Sapphire lasers 
transfer-cavity locked to a Rb-stabilized diode laser. 

As the 421-nm laser system is detuned to the red of 
isotopes 160 through 164, steady state MOT popula- 



tions form in proportion to natural abundance, except 
for the fermions ^^^Dy and ^^^Dy. These isotopes' MOT 
populations are 5/6 x and 1/6 x the expected, respec- 
tively, which is likely due to poor optical pumping to the 
F = 21/2 state [IS], where J, J, and F = J / are the 
total electronic, nuclear, and total angular momenta. 

Figure [2][ a) describes our model of the Dy MOT recy- 
cling mechanism, which is a refinement of that proposed 
for the Er MOT [IB]. The blue cooling laser excites a 
fraction /ox < 1/2 of the population to the 421-nm level, 
where it can decay with branching ratio B = Ri/T to 
the metastable states. Upon decay, a fraction pf^x of 
the population is captured in the magnetic quadrupole 
trap (MT) of the MOT at rate pfe^Ri- MOT recycUng 
data, examples of which are shown in Figs.|2][b) and (c), 
support a dual decay path through the dark metastable 
states; with probability q {1 — q) the population decays 
through the slow (fast) branch at rate i?siow (-Rfast)- Pop- 
ulation that reaches the ground state reloads the MOT 
at rate i?rcioadi which depends on MOT parameters. 

There are two loss rates from the otherwise closed 
system (once loading from the Zeeman slower ceases). 
A portion of the MOT population can decay to 
non-magnetically trapped metastable states at rate 
_RiossMOT = (1 — p)fcxRi- Additionally, metastable and 
ground state MT populations can be lost due to back- 
ground and two-body inelastic collisions. Two-body loss 
is difficult to quantify in the metastable states, but we 
investigate ground state MT loss below: i?siow is faster 
than the ground state component of i?iossMT , and we ne- 
glect -RiossMT in the rate equations below. 

Following the procedure outlined in Ref . ^16j , the rates 
in the MOT recycling model are determined by fitting 
the following equations to sets of MOT decay transients 
taken at various combinations of MOT beam intensity 
and detuning (see as an example the data in Fig. [3f[b)): 

A^MOT = ^reload A^MT " /cx-Rl A'^MOT, 

A'fast = (1 - q)p.f ckRiNmOT - -RfastA^fast, 

A'slow = qpf cxRiNmOT ~ ^slowA"slow, 

AmT — ^fastA^fast + -RslowA"slow — -R-rcload A^MT , (1) 

where Amot, A^fast, A^siow, and A^mt are the popula- 
tions of the MOT, fast (slow) metastable state decay 
channel, and MT, respectively. To obtain MOT recy- 
cling data with well defined initial conditions, the MOT 
is loaded with the Zeeman slowed atomic beam until 
reaching steady state, then the MOT, slower, and atomic 
beams are extinguished for = 1 s, during which most 
of the population in the metastable states decay to the 
ground state. The population equations are numerically 
fit to the data for nine combinations of MOT power and 
detuning for the fermion ^^'^Dy and eight for the boson 
^^■^Dy. We verified that the hold time tm in the MT is suf- 
ficiently long that any residual population in Ngiow does 
not affect fit results; all other variables are left free to 
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vary. Ri is extracted from the product R = fcxRi by si- 
multaneously fitting the R's for each isotope to the func- 
tion R{I,S) = i?i//(2 + 21 + 2{2S/T)^), where / = I/I^ 
and Is « 2.7 x 58 mW/cm^ in the MOT. 

Averaging the results, we find that for both isotopes 

[rI'^RZIRT^] = [1170(20), 19(2), 1.5(1)] s-\ 
[i?r,^fa't,^s\ot] = [1700(100), 29(1), 2.3(1)] s-\ 

and [p,q] = [0.82(1), 0.73(1)]. The i?}64 ^.^^^g jg consis- 
tent with the corresponding bosonic ^^^Er (/ = 0) MOT 
quantity, and we suggest that the isotope-induced differ- 
ence in Dy MOT decay rates arises from ^^'^Dy's hyper- 
fine structure. Dy Zeeman slowing and MOT collection 
is possible because of the small B ^ 7 x 10~^ [21 . The p 
and q values indicate that 82% of the atoms are captured 
by the MT and 73% of those cascade through the slow 
channel. Dy level linewidths are mostly unknown, but 
we speculate that much of this population is captured 
by the low-lying and long-lived states near the 1322-nm 
telecom transition. 

Using the measured i?siow rates as fixed constants, data 
of the maximum MOT recapture population (peaks of 
data such as in Fig. |2jb)) versus tm are fit to an expo- 
nential with only the initial and final MOT population as 
free variables. Good fits are obtained for < 1 s data, 
as shown in Fig. ^c); we plot the fitted curve out to 
later times for comparison to data [52] . From these plots 
we can see that the population in the steady state MOT 
is only a fraction of the total number of trapped atoms. 
The "hidden" population in the continuously loaded MT 
is several times larger, filling the ground state MT at 
rate Rsiow With respect to the steady state MOT, we 
measure nearly 2.5 x (3.5 x) more ^^^Dy (^^^^Dy) atoms 
in the recaptured MOT. At optimal MOT parameters 
and after = 2 s of MT loading, we measure the to- 
tal population of laser cooled and trapped atoms to be 
5 X 10^ for i64Dy and 3 x 10^ for ^^^Dy. Populations are 
measured via fluorescence collection on a fast photodiode 
with atom number calibration to < 10% from absorption 
imaging [12]. 

Population loss at times t,n > 1 s is likely due to back- 
ground and inelastic collisions of ground state atoms in 
the MT. (Majorana spin flip losses are negligible at these 
tm's and trap densities.) We investigate these collision 
rates by absorption imaging the atoms confined in the 
MT at various delay times. After loading the MOT, 
we extinguish the MOT, Zeeman slower, and atomic 
beams while maintaining a constant Vi3. The MT is 
turned off at tm, and after 1 ms, an absorption image 
integrates on a 16-bit CCD camera for 200 /iS. Account- 
ing for gravity —gz, the images taken in the x-z plane 
containing the quadrupole axis (along z) are fit to the 
functional form |23j : 

p,(x) cx exp(-2|a:|yr^/x)(l -h 2|x|yr^/x), 
p,{z) cx exp(-2|z|/z - 2gz/z){l + 2\z\/z), 




Time 1^ in magnetic trap (s) 

FIG. 3. (color online), (a) Decay of magnetic trap mean 
density for three gradients. Density profiles are consistent 
with a single thermal distribution after tm = 2 s. Data are fit 
to obtain one (7) and two-body (/3) collisional loss rates, (b) 
Heating in VzB = 22.2 G/cm magnetic trap, (c) Increase in 
mean magnetic moment p, for the same trap. An exponential 
fit gives a fl{t) time constant of r = 0.22(5) s~^. 

where g = Mg / {pSI yB), ft is the mean magnetic mo- 
ment of the atomic cloud, and x = 2kBT / {fiS/ xB) (and 
similarly for z). We extract the mean cloud density n, 
temperature T, and Ji versus t,„, which we plot in Fig. [3] 
for several MT gradients. 

The -'^^''^Dy density decay data in Fig.[3]|^a) are fit to the 
one plus two-body decay equation n=—jn—/3n^, where 
7 is the loss rate due to background collisions and /3 is 
a measure of inelastic two-body losses. For the three 
trap gradients, the extracted /3's and 7's are consistent 
with one another within la; the weighted means are 
[/3,7] = [2.1(2)xl0-i" cmVs, 5.6(3) x 10-2 s~^]. Fitting 
the data to n = —fSn^ results in a worse ■ Temperature 
heating data — e.g.. Fig. |3jb) — are fit to an exponential, 
resulting at early times to a heating rate of 17(2) fiK/s, 
which is consistent with the heating rate in spin un- 
polarized Er and Cr MTs [H [H] and is likely due to 
spin-relaxation collisions. The MT population is initially 
distributed among the weak-field seeking Zeeman states, 
and spin exchange collisions tend to polarize the sam- 
ple toward mj = 8. Indeed, fl increases with time, see 
Fig. [3]Jc), and fi{t) reaches 8/15 within t„j = 4 s in the 
high gradient VzB = 34.7 G/cm trap. 

The measured (3 is likely due to an unresolved combi- 
nation of inelastic spin exchange, magnetic dipole-dipole 
relaxation (MDDR), and anisotropic electrostatic-driven 
spin relaxation collisions. Compared to (3ci ~ 3x10^^^ 
cm'^/s in 200 /iK spin-polarized ^^Cr [53|, (3 is consis- 
tent with a MDDR scaling c>c n'^, though several times 
smaller than the full inelastic MDDR scaling presented in 
Ref. [2?. In addition, (3 is consistent with non-maximally 
spin-polarized Cr (LlxlO"^" cm^/s) [24j and with 500 
mK inelastic Er-Er and Tm-Tm collision rates — 3.0 and 
1.1x10-^" cm^/s, respectively [T5]; the large magni- 
tude of the latter rates are attributed to anisotropic 
electrostatic-driven spin relaxation collisions of these lan- 
thanides. Further measurements will aim to elucidate 
these collisional processes in this complex atom [TO] . 

Temperature and density profiles contain unusual fea- 
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FIG. 4. (color online), (a)-(c) ^''^Dy absorption images at 
time td after release from MOT revealing a dense stripe with 
anisotropic temperature, (d)-(f) Column density integrations 
along z versus x. Fits to the data (blue) are to a double 
(white) and single (red) Gaussian. Additional images in the 
x-y plane reveal the stripe's azimuthal symmetry. 

tures common to both ^^'*Dy and ^^^Dy MOTs. When 
care is taken to retroreflect and power balance all pairs 
of MOT beams, the MOTs are typically 30% fewer in 
population, though have larger mean density (10^^ ver- 
sus 10"'^° cm"'^) than those with slight misadjustment. In 
ballistic expansion after extinguishing the MOT, a dual 
component gas is observed comprised of a dense symmet- 
ric core of ~200 /iK atoms surrounded by a hot, 2-3 mK 
shell containing 70% of the atoms. Doppler coohng the- 
ory in ID predicts a cloud temperature of 1.2 mK for the 
MOT parameters, but comparisons to the Er MOT [23j 
suggest the entire intra-MOT population should be sub- 
Doppler cooled to ~100 fxK as a consequence of the near 
equal Lande (/-factors in the ground and excited states 
{Ag/g = 1.7%). Despite Ag being slightly larger than in 
Er, numerical ID sub-Doppler cooling simulations (based 
on Ref. [23j and references within) indicate that the entire 
Dy MOT — accounting for typical size and VB — should 
be sub-Doppler cooled. The Er MOT, limited by Zeeman 
and MOT laser power, contained 500 x fewer atoms; it is 
possible that the hot shell forms subsequent to the cold 
inner core, which will be investigated. At t,n — 2 s, MT 
temperatures are 100-300 /iK depending on VS. 

While ultracold cores are observed in other MOTs [2S] , 
a unique sub-Doppler cooled structure forms when the 
Dy MOT beams are slightly misbalanced or misaligned; 
an example of which is shown in Fig. [4] Depending on the 
particular misadjustment, either a dense vertical or hor- 
izontal stripe (or both) appear in the core of the cloud. 
The outer hot atoms remain at = 2-3 mK temper- 
atures, but the (vertical) stripe population acquires an 
anisotropic temperature distribution of T^' « and 
< 10 /iK. This latter temperature is consistent the ID 
numerical sub-Doppler cooling calculation, but its precise 
measurement is hampered by the hot cloud presence. We 
offer no explanation other than to note no detection of 
stripe spin-polarization in Stern-Gerlach measurements 
(again hampered by the rapid hot cloud expansion), and 
that rf T" as VB decreases [19]. Under optimal 



operating conditions, the Dy MOT phase space density 
is as large as 10"''. 

Future work includes narrow- line cooling to ^^1 /iK on 
the 1001-nm line [25] and loading into a crossed opti- 
cal dipole trap. Once optically confined, the large MT 
losses measured here might be avoided by trapping in 
the lowest energy Zeeman state. Elastic and inelastic 
collision rates in the single partial wave regime will be 
measured before attempting to evaporatively cool to de- 
generacy. The ultracold Dy produced with this method 
will open new avenues for research in quantum gases, pre- 
cision measurement, and quantum information science. 
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